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ABSTRACT 

When we model black hole accretion sources such as active galactic nuclei 
and black hole X-ray binaries as advection-dominated accretion flows (ADAFs), 
it is neccesary to use the global solution to the equations rather than the simpler 
self-similar solution, since the latter is inaccurate in the region near the black 
hole where most of the radiation is emitted. However, technically, it is a difficult 
task to calculate the global solution because of the transonic nature of the flow, 
which makes it a two-point boundary value problem. In this paper we propose 
a simplifled approach for calculating the global ADAF solution. We replace the 
radial momentum equation by a simple algebraic relation between the angular 
velocity of the gas and the Keplerian angular velocity, while keeping all other 
equations unchanged. It is then easy to solve the differential energy equations 
to obtain an approximate global solution. By adjusting the free parameters, we 
find that for almost any accretion rate and for a = 0.1 — 0.3 we can get good 
simplified global solutions. The predicted spectra from the approximate solutions 
are very close to the spectra obtained from the true global solutions. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 
active — galaxies: nuclei — hydrodynamics 

1. Introduction 

Advection-dominated accretion fiow (ADAF) is an important type of solution for black 
hole accretion. A prominent feature of an ADAF compared to the standard thin disk is its 
low radiative efficient at low accretion rates (Narayan & Yi 1994, hereafter NY94; Narayan 
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& Yi 1995; see Narayan, Mahadevan & Quataert 1998 and Kato, Pukue & Mineshige 1998 
for reviews). 

The ADAF solution has received much attention in the past years because it successfully 
explains why some nearby galaxies are so dim even though their accretion rates are not very 
small (see Narayan 2005, Yuan 2007, and Ho 2008 for reviews). The best evidence comes 
from the supermassivc black hole in our Galactic center, Sgr A* (Narayan, Mahadevan & 
Yi 1995; Manmoto, Kusimoze & Mineshige 1997; Yuan, Quataert & Narayan 2003). From 
Chandra observations and Bondi theory, we can estimate the mass accretion rate of Sgr A*. 
If the accretion flow were not an ADAF but a standard thin disk, the luminosity would be 
five orders of magnitude larger than observed (Yuan, Quataert & Narayan 2003). 

Another attractive feature of an ADAF is that it can partly solve the problem of the 
origin of X-ray emission from accretion flows. The temperature of a standard thin disk 
at the inner disk is only ~ lO^K for a supermassivc black hole, too low to produce X-ray 
emission (Frank, King & Raine 2002). A hot corona has been thought to be responsible 
for the X-ray emission, but recent MHD simulations of disks show that they have hardly 
any coronae (Hirose, Krolik & Stone 2006). On the other hand, the temperature of an 
ADAF is high enough to produce X-ray emission. Of course, a canonical ADAF exists only 
below a critical accretion rate Afcrit ~ a^^Edd with MEdd = lO^Edd/c^ and a is the viscous 
parameter, which corresponds to ~ (3 — 4)%LEdd at most (Esin, McClintock & Narayan 
1997). Therefore it cannot explain luminous X-ray sources such as quasars. The luminous 
hot accretion flow (LHAF; Yuan 2001), which is an extension of an ADAF to higher accreton 
rates, is promising, but the details of this model have not been fully worked out (see Yuan 
et al. 2007 for an example of application to luminous black hole X-ray binaries). 

In spite of the great success of ADAFs, more work is required to test the model. On the 
one hand, it would useful to expand the application of ADAFs to more sources, and on the 
another hand, such modeling is expected to help us understand some important microphysical 
issues which are still unclear. One example of the latter is the value of 5 (defined in eqs. 4 and 
5), which measures the amount of direct electron heating through viscous dissipation in a hot 
accretion fiow (Sharma et al. 2007). Another is the potential importance of collective plasma 
effects which will determine how realistic the two-temperature assumption is (Begelman & 
Chiueh 1988; Yuan et al. 2006). 

For such work, the global solution rather than the self-similar solution of the ADAF 
equations is required. This is because most of the radiation of an ADAF comes from its 

innermost region where the self-similar solution breaks down. However, it is technically very 
difficult to calculate the global solution of an ADAF. An ADAF is transonic, and thus its 
global solution should satisfy the sonic-point condition in addition to the outer boundary 



-3- 



condition. Mathematically, it is a two point boundary value problem and not easy to deal 
with. This is an obstacle to the wide application of the ADAF model. 

In this paper we propose a simplified global ADAF solution. We adopt a simple algebraic 
relation to replace the radial momentum equation, thus avoiding the two point boundary 
value problem. We present our approach in §2 and show some examples in §3. The final 
section is devoted to a short summary. Watarai (2007) has recently presented related work. 



2. The simplified global ADAF model 

The basic equations of an ADAF describe the conservation of mass, radial and aximuthal 
components of the momentum, and energy (e.g., Narayan, Mahadevan & Quataert 1998): 



M = -AnrHpv (1) 

,- = in^-al)r--- (2) 

v(Qr^ — j) = —arcl (3) 

pvT,-j^ = (1 - S)q+ - (4) 

dSe 

PvTe-T^ = qie-q - (5) 



All the quantities have their usual meaning. In the present paper we do not include outfiows 
from the ADAF, but it is easy to extend our calculation to that case by simply using a radius- 
dependent mass accretion rate, M = Mo(r/rout)* with s > being a constant (e.g., Yuan, 
Quataert & Narayan 2003). The quantity 6 in equations (4) and (5) describes the fraction 
of the turbulent dissipation rate which directly heats electrons; we set 6 = 0.3. The 
quantity Qie describes the energy transfer rate from ions to electrons by Coulomb collision, 
and q~ is the radiative cooling rate. We consider synchrotron and bremsstrahlung emissions 
and their Comptonization. The details of the calculation of the spectrum can be found in 
Yuan, Quataert & Narayan (2003). We consider a Schwarzschild black hole and adopt the 
Paczyhski & Wiita (1980) potential to mimic its geometry. 

The most difficult part of solving the global solution is the radial momentum equation 
(2). Our key idea of simphfying the global ADAF solution is to replace this differential 
equation by the following simple algebraic relation: 



(6) 
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with 

/o for r > rms 

where rms — Sr^ = 6GM/c^ is the innermost stable circular orbit. 



r < r„ 



The above simplification is based on the following physical consideration. The imme- 
diate idea we think of to simplify the radial momentum equation is to use the self-similar 
solution obtained by NY94. Consider eqs. (7)-(9) in NY94. We can use eq. (9) to solve for 
e' in terms of the sound speed c^: 

cj cj 2 

Wc can then substitute this in eqs. (8) and (7) in NY94 to obtain for a fully advection- 
dominated flow (e' = e) 

Q / bclV'^ 10-67 , 

7^ = 1 - TT-T- = -T. 7- = const. (8) 

V ^^l) 97-5 ^ ' 

Here 7 is the adiabatic index. We therefore in principle could set / in eq. (6) to this constant. 
However, we find that the simplfied solution is very sensitive to the value of /. The reason 
is that, as we will see, the radial velocity is sensitive to the value of / (ref. eq. 9 below). 
The velocity determines the density, and also the temperature via the energy equations, two 
quantities that determine the emitted spectrum. We therefore set / as a free parameter 
which we adjust for different accretion parameters M and a to get the best approximation. 

Because the angular momentum Qr^ in a global solution keeps decreasing with decreas- 
ing radius, while the Keplerian angular momentum Q^xf^ begins to increase when r < r^g 
(Fig. 1), / cannot be a constant when r < Tms- Instead we require the angular momentum 
to be continuous at r^s and assume that it is proportional to {r /rms)"'- After some tests we 
set n = 0.5, independent of the values of M and a. Thus n is not a free parameter in our 
model. 

Substituting eq.(6) into eq.(3) we have 

(9) 



The quantity j is the specific angular momentum of the accretion gas when it falls into the 
black hole and it is the eigenvalue of the exact global solution. In our simplified model, we 
set j as the second free parameter and adjust its value to get the best approximate solution 
for V. 
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Substituting eqs. (1), (6) and (9) into the energy equations for ions and electrons, eqs. 
(4) and (5), we have two differential equations with two unknown variables, Tj and Tg. All 
other quantites such as v, p, and H can be expressed as simple functions of Tj and T^. for 
a given M and a and assumed values of the free parameters /o and j. When Tj and Tg are 
given at the outer boundary, we can easily integrate the differential equations inwards to get 
the approximate global solution. 

We adjust the values of /o and j for different M and a to get the best simplified global 
solution. Here "best" means that the profiles of all quantities such as p, f , Tg, Tj, and most 
importantly, the emitted spectrum, are very close to the exact global solution. Because our 
main purpose is to model the continuum spectrum of black hole sources (AGNs and black 
hole X-ray binaries), our first priority will be the closeness of the spectrum when we judge 
how good a simplified solution is. 

In the calculations presented here, we set M of the approximate solution equal to M 
of the global solution it is meant to fit. However, in real applications, we only know the 
spectrum rather than M. So it might be more reahstic to adjust Mguap of the simplified 
solution to fit the spectrum produced by the exact global solution with a given Mexact rather 
than setting Mgimp = -/^^exact- Fortunately we find that Mgimp and Mexact are very close, 
typically eviating by no more than ~ 3%. 

3. Results 

When modehng black hole sources with an ADAF model, the accretion rate M spans 
a wide range, say from 10~^MEdd to 10~^MEdd- But the value of a adopted in ADAF 
modehng (e.g., Narayan, Mahadevan & Quataert 1998) is usually within a very narrow 
range, a = 0.1 — 0.3. This is also supported by MHD numerical simulations of accretion 
flows (e.g., Hawley & Krolik 2001). We adjust the values of j and /q to obtain the "best" 
simplified ADAF solutions for M and a within the above ranges. As state below we find 
that the same set of (j, /o) often holds for quite a wide range of M. 

3.1. a = 0.3: /o = 0.33, j = 0.98 for any M 

We first present results for a — 0.3. We find that in this case the simplified global 
solution with /o = 0.33, j = 0.98 gives a satisfactory spectrum for any M. Figs. 1 & 2 
give two examples with M = 10^^ and lO^^Mgcid) respectively. The dashed lines in the 
figure denote the exact global solution while the solid lines are for the simplified global 
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solution. The plots in each figure show the emitted spectrum, Mach number, electron and 
ion temperature, density, and the angular momentum, respectively. For Fig. 1, the outer 
boundary is at lOV^ and the outer boundary condition is Tj = 0.2Tvir, Tg = 0.19Tvii. with 
the virial temperature T^ij- = 3.6 x 10^^ (r^/r). For Fig. 2, the outer boundary is at lO^r^ 
and the outer boundary condition is Tj = O.GTyir, Tg = O.OSTyir. 

3.2. a = 0.1: /o = 0.33, j = 1.08 for M < lO-^MEdd 

When a = 0.1, it is hard for a single set of (j, /o) to give a good solution for all M. 
When M is relatively low, M < 10~^MEdd, we find /o = 0.33, j = 1.08 gives a satisfactory 
solution. Fig. 3 shows an example with M = 10~^MEdd- The outer boundary conditions are 
Ti = 0.2rvir, Te = O.lQTvir at rout = 10%. 

3.3. a = 0.1: /o = 0.15, j = 0.49 for M > 5 x lO-^Medd 

When M is relatively high, M > 5 x 10~^MEdd, we find /o = 0.15, j = 0.49 gives a 
satisfactory solution. Fig. 4 shows an example of M = lO^^Mgdd- The outer boundary 
conditions are Tj = O.GTvu-, Tg = 0.08Tvir at rout = lO^r^. Note this solution is in the 
regime of LHAF, because 10~^MEdd is well above the critical accretion rate of an ADAF 

Merit « lO-^MEdd- 



3.4. Other values of a 

For other values of a, we find that simply using the "linear fit" values of (j, /o) between 
those for a = 0.3 and 0.1 gives a good solution. For example, for a = 0.2 and M = lO^'^MEdd, 
the values of j and /o are just /o = (0.33 + 0.33)/2 = 0.33, j = (0.98 + 1.08)/2 = 1.03. For 
a = 0.2 and M = 0.08MEdd, /o = (0.33 + 0.15)/2 = 0.24, j = (0.98 + 0.49)/2 = 0.735. 

4. Summciry 

The global solution of ADAFs is difficult to calculate because it is mathematically a 
two point boundary value problem. This hampers wide application of the ADAF model. We 
propose a simplifed global solution to overcome this difficulty. Prompted by the self-similar 
solution of ADAFs, we replace the radial momentum equation, which is the most difficult to 
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handle, with a simple algebraic relation (eq. 6) , and then solve the remaining two diferential 
equations (eqs. 4 & 5). We adjust the two free parameters (j and /o in eqs. 6 and 9) so 
that we obtain the best approximation compared to the exact global ADAF solution. The 
spectra of the simplified solutions are impressively good, as shown in Figs. 1-4 for various 
M and a. 

We have been unable to identify a single set of values of the two adjustable constants 
j and /o which works for all ADAF models. However, when a is large, say ~ 0.3, we find 
that j = 0.98, /o = 0.33 gives very good results for all accretion rates M for which an 
ADAF solution is possible. Recent work by Sharma ct al. (2006) suggests that the viscosity 
parameter in the coUisionlcss plasma in a hot accretion flow will be larger than in a standard 
thin disk. Therefore, a ~ 0.3 is probably not unrealistic for an ADAF. It would thus be 
reasonable to use a single set of parameters, a — 0.3, j — 0.98, /o = 0.33, for practical 
applications of the approximate global model described here. 

The success of the present work encourages us to extend our approach to the case of 
a slim disk, which is an extension of the standard thin disk to accretion rates above the 
Eddington rate (Abramowicz et al. 1988). It potentially has important apphcation in ULXs 
and narrow-hne Seyfert 1 galaxies (Mineshige et al. 2000; Watarai et al. 2001). We hope to 
report the results in a future paper. 

This work was supported in part by the Natural Science Foundation of China (grant 
10773024), Shanghai Pujiang Program, Baircn Program of CAS (F.Y.), the Knowledge In- 
novation Program of CAS, and Shanghai Postdoctoral Scientific Program (R.M.). 
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Fig. 1. — The simplified (solid) and exact (dashed) global ADAF solutions for M 
10~^MEdd) = 0.3. The parameters are /o = 0.33, j = 0.98. 




Fig. 2. — The simplified (solid) and exact (dashed) global ADAF solutions for M 
10~^MEdd, CK = 0.3. The parameters are /o = 0.33, j — 0.98. 
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Fig. 3. — The simplified (solid) and exact (dashed) global ADAF solutions for M — 
10~^MEdd, CK = 0.1. The parameters are /o = 0.33, j — 1.08. 




Fig. 4. — The simplified (solid) and exact (dashed) global ADAF solutions for M — 
10"^MEdd, a = 0.1. The parameters are /o = 0.15, j = 0.49. 



